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Abstract 

Finite temperature density functional theory requires representations for the internal energy, 
entropy, and free energy as functionals of the local density field. A central formal difficulty for 
an orbital-free representation is construction of the corresponding functionals for non-interacting 
particles in an arbitrary external potential. That problem is posed here in the context of the 
equilibrium statistical mechanics of an inhomogeneous system. The density functionals are defined 
and shown to be equal to the extremal state for a functional of the reduced one-particle statistical 
operators. Convexity of the latter functionals implies a class of general inequalities. First, it is 
shown that the familiar von Weizsacker lower bound for zero temperature functionals applies at 
finite temperature as well. An upper bound is obtained in terms of a single-particle statistical 
operator corresponding to the Thomas-Fermi approximation. Next, the behavior of the density 
functionals under coordinate scaling is obtained. The inequalities are exploited to obtain a class of 
upper and lower bounds at constant temperature, and a complementary class at constant density. 
The utility of such constraints and their relationship to corresponding results at zero temperature 
are discussed. 
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I. INTRODUCTION 



Density functional theory (DFT) is an important structure for posing problems of equilib- 
rium many-body physics in a non-perturbative context. Most applications, with many signif- 
icant successes, have been in the extremes of zero-temperature ground state systems [l-7| or 



high-temperature classical systems jSHlOj. Shortly after the ground-state Hohenberg-Kohn 



Sham formulation, a finite-temperature generalization was given [ll|] but it has attracted 
far less development than the ground-state version 12|. Recently, attention has been di- 
rected toward systems in a difficult intermediate domain known as "warm dense matter". 
Extensions of existing density functional approaches (including procedures to develop ap- 
proximate functionals non-empirically) to this domain have met with challenges. One is 
the temperature-dependence itself of the functionals. Another is technical, associated with 
solving a self-consistent single-particle eigenvalue problem (the Kohn-Sham equation) for 
which a growing number of eigenvalues and eigenf unctions (orbitals) is required as the tem- 
perature is increased at solid densities. For many applications, such as determining forces 
for Born-Oppenheimer molecular dynamics of ions in warm, dense matter, this growth pro- 
duces a computational bottleneck. The Kohn-Sham procedure is not intrinsic to DFT, so an 
alternative approach is to return to the expression of the entire free energy as a functional 
of the density. Once known, the extremum condition of DFT becomes a single non-linear, 
non-local equation for the density whose solution replaces that constructed from the Kohn- 
Sham orbitals and eigenfunctions. Solving that Euler equation is facilitated by introducing 
the Kohn-Sham separation into a non-interacting functional, which is presumably dominant, 
and the remainder, but without introducing orbital-dependent expressions. This approach 
is known as orbital free DFT (OFDFT). An example is given by the Thomas- Fermi-Dirac 
model, which is a local density approximation. 

A critical component for OFDFT is the free energy density functional for an inhomoge- 
neous but non-interacting system. The objective here is to define the free energy, energy, 
and entropy density functionals for a non-interacting system at finite temperature in forms 
suitable for constructing approximations. (Note the correspondence with the ground state 
theory, to wit, the definition of the Kohn-Sham kinetic energy and exact exchange energy.) 
To aid in that construction important exact constraints are obtained here as well. These 
include scaling laws, upper and lower bounds, and inequalities. Development and use of such 
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constraints is promising in light of the history of analogous developments in zero-temperature 
DFT. 

The approach here is to exploit the relationship of DFT to the equilibrium statistical 
mechanics of non-uniform systems 3]. The variational context of DFT is simply related 
to the extremal properties of the ensembles representing equilibrium systems. In the next 
section, thermodynamics is described in the grand ensemble, for which the thermodynamic 
variables are the temperature T and chemical potential /i. The presence of an external single- 
particle potential v{r) transforms the chemical potential to a local function /i (r) = ^ — v{r). 
A change of thermodynamic variables T, fi (r) to temperature and density T, n (r) is effected 
by an appropriate Legendre transformation. The resulting thermodynamic functionals of 
T, n (r) are the desired density functionals for DFT. The extremum condition of DFT is 
simply the thermodynamic identity relating the functional derivative of the free energy to 
the chemical potential. 

This context for DFT is described briefly in the next section and then specialized to the 
non-interacting system in section lllli The equation for the density provides the definition 
of the local chemical potential as a functional of the density. The free energy, energy, and 
entropy are represented in terms of the equilibrium single-particle reduced density matrix, 
or Fermi operator. The latter depends on the chemical potential, so that in this way these 
properties are identified as functionals of the density through the chemical potential. Next, 
corresponding functionals of an arbitrary single-particle statistical operator are defined. 
Using convexity of the free energy functional of statistical operators, it is shown that the 
minimum occurs for the Fermi operator, resulting in the free energy density functional being 
a minimum for the free energy statistical operator functionals. This is the source for the 
inequahties described subsequently. 

In section IIVI the transformation of the density functionals under coordinate scaling is 
obtained from dimensional analysis. An alternative constructive derivation which connects 
explicitly with ground-state uniform scaling 15l-ll7l] is given in Appendix [Bl Next, the 
familiar von Weizsacker lower bound on the kinetic energy for zero-temperature DFT 

0- 



20[ | is shown to hold for finite temperature as well. A statistical operator underlying the 



Thomas-Fermi approximation for the density is defined via an appropriate Wigner function 
representation and used to obtain an upper bound for the free energy. Finally, two classes of 
inequalities are established for all density functionals using the scaling laws. One is for fixed 
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temperature at different densities, while the other is for fixed density and different temper- 
atures. For the case of differential scaling, the results reduce to differentials of equilibrium 
thermodynamics. 

In the last section these results are discussed and their utility for constraint-based ap- 
proximations described. 



II. EQUILIBRIUM STATISTICAL MECHANICS FOR A NON-UNIFORM SYS- 
TEM 



While treatment of a homogeneous non-interacting Fermion system is a textbook staple, 
the statistical mechanics of m-homogeneous non-interacting systems is less well documented. 
The formal structure of ensemble DFT has been discussed in diverse forms 1211 25 1 since 
Ref. 111. Some of those formulations resemble the ground-state theory more than others. 
Homogeneous scaling at finite temperature has been treated recently in a manner closely 
related to ground-state scaling 26|. The present approach is complementary, in the sense 
of providing a close, systematic connection with statistical mechanical methods widely used 
in the treatment of electronically hot systems, e.g. plasmas. To set notation and make the 
context explicit, we begin with a brief summary of the relevant statistical mechanics. 

The system of interest is comprised of electrons in a uniform, rigid, neutralizing back- 
ground ("one component plasma", "jellium", "electron gas") I27I], and in an external poten- 
tial V. The Hamiltonian is 

N 

H = He + J2<^o.). (1) 



A caret over a symbol indicates that it is an operator. Here 

H - V^ + i V - V /rfr^^ + i ! dvdv'^^ (2) 

^<^-Z.2m„ + 2 lq.-q,| ^^|q„-r|+^y ^"^"|r-rr 

The integrations extend over the volume V and nf, is the average number density of the 
background chosen for over-all charge neutrality. The first two terms are the kinetic energy 
and electron-electron Coulomb interaction energies. The third term is the interaction energy 
of the electrons with the neutralizing background (charge density rifoe), while the last term is 
the interaction energy among elements of that uniform background. The origin and necessity 
for this neutralizing background arises from the long range nature of the Coulomb potential - 
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a macroscopic equilibrium state for the electrons alone does not exist [28!]. Charge neutrality 
is a necessity for thermodynamic stability, e.g. systems of both electrons and positive ions. 
The system considered here is a model for electrons in such two-component systems, where 



the uniform, rigid, neutralizing background assures that thermodynamics exists 
The electron density operator is defined as 

N 



29 



30|. 



[r) = ^5(r-q„). (3) 

Note that 



a=l 



TV 



n(r)n(r') = ^ 5 (r - q„) (r' - %) + 5(r - r')n(r) (4) 
Then Eq. ([T]) becomes 

N ^2 r 

^2m, J 
r 

+ I y rfrrfr'^— [(n{r) - rib) {n{r') - rib) - 5(r - r>b] (5) 

The problem of interest is to determine the equilibrium electron density n{r) and the ther- 
modynamics for this inhomogeneous electron gas. 

The thermodynamic properties for this class of inhomogeneous systems are determined 
in the grand canonical ensemble as functions of the inverse temperature /3 = l/ksT and 
functionals of the local chemical potential 

fx{r) = fx-v(r). (6) 

which combines the external potential v{r) with the constraint on average total particle 
number, the constant chemical potential /i. The fundamental thermodynamic potential for 
this ensemble, fl{f3 \ /i), is related simply to the pressure 

oo 

Q{I3 I ^) = -p{(3 I ^i)V = In J2 Tr(^)e-^(^~^*'^('")"('")) . (7) 

N=0 

The trace extends over the Hilbert space for electrons. The notation {P \ fi) indicates 
variable dependence and functional dependence, respectively. The thermodynamical average 
electron density as a function of (3 and functional of /i(r) follows directly as 

nM|,) = -^^ . (8) 
d/i (r) /3 



Similarly, the internal energy is determined from the derivative with respect to f3 

Now consider a change of variables from I3,fi{r) to f3,n{r) by a Legendre transformation 
to the free energy 

F(,|,,,)^.(,|,) + /*,W„W. (10) 

By construction F{P \ n) is independent of /i(r) {i.e., SF{P \ n)/6fi{r) = 0). Therefore it is 
"universal" in the sense that the functional F{P | ■) is the same for all external potentials, 
with its argument evaluated at the corresponding different density fields. In all of the 
following, only density fields n(r,/3 | fi) that can be represented by a chemical potential 
yu(r) as in ([7]) and ([8]) are considered. These are the /i-representable physical densities of 
thermodynamics. The corresponding ground-state theory requirement on the density is v- 
represent ability. With this restriction on the class of densities, Eqs. ([7]) and (jH]) provide the 
constructive definition for the free energy functional: calculate Q{(3 \ fi) as a functional of 
H from its explicit definition; invert Eq. ([8]) for the chemical potential as a functional of the 
density fi{r, (3 \ n); substitute fi{r, f3 \ n) into the right side of ( fTOj) . 
The functional derivative oi F{(5 \n) gives the inverse of ([8]) 

5F{P I n) 

n) = (11) 

on (r) 13 

This is an identity which expresses the one - one relationship of a /i (r) with a unique density 
n (r). However, it also can be viewed as an equation to find the density associated with a 
given fi (r). For example, consider the case where /io (r) is given, and the objective is to find 
the corresponding uq (r) . This can be accomplished by solving the equation 

6F{/3 I n) 



(12) 



Sn (r) 

The unique solution is n (r) = uq (r). This is no longer an identity, but an equation to be 
solved: equality holds only for the specific density uq (r). 

Up to here the development is simply equilibrium quantum statistical mechanics. But the 
same result follows from DFT in an equivalent way. In terms of the stipulated local chemical 



potential /io(r), ensemble constrained search [3l| enables us to define the functional 

{P\n) = F{/3 \n)- j drfio (r) n(r). (13) 



The condition for an extremum of {f3 \ n) gives (fT2l) as its Euler equation 

Sn^,{f3\n) _^ _ , ._6F{f3\n 



From Eqs. note that 

n^,{/3 \n) = n{/3\fL)- J dr (/io (r) - fi{r,/3 \ n)) n{r), (15) 

and so 

n^,{(3\no) = n{f3\fio). (16) 
The case of primary interest is an external potential due to a given configuration of ions 

^;o(r) = /i - /io(r) = Ze' J eirdr^ ^^^"^^ n,(r) = f; 5 (r - R..) (17) 



(14) 



a=l 



where Z is the ion charge number, and ZNi = N. For simplicity, we discuss only a single 
ionic species; generalization is obvious. 

Determination of the electron density field for given ion configurations (hence, given /io(r)) 
involves constructing F(/3 | n) and solving ( 1T4|) for no(r). The thermodynamic properties 
then are determined from F{P \ rio). The problem of constructing F{P \ n) is separated 
into two parts: finding the corresponding functional for non-interacting electrons (but in 
the external fields of the ions), F^^\(3 \ n), and finding the functional associated with the 
contribution from Coulomb interactions, AF(/3 | n). This is a logical separation, since 
F^^\(3 I n) is itself a universal functional, with properties similar to F((3 \ n). For example, 
there is a unique functional fj.^^\r,(3 \ n) associated with F^^\l3 \ n) for each n given by 



5F(o)(/3 



n] 



on (r) 13 

The notation /x^*^^ specifies that the functional is that for the non-interacting relationship, 
with the same chemical potential and external potential 

/i(o)(r,/3 I n) = /i-t;(r) . (19) 

The Coulomb contributions AF(/3 | n) ("excess free energy" in statistical mechanical 
usage) conventionally are divided into a mean-field (Hartree) contribution plus a remainder, 
the exchange-correlation contribution 

AF(/3 \n) = \j dvdr'n (r) ^^^^^n (r') + F^,{I3 \ n), (20) 



so that the total free energy is composed as 

F{f3 I n) = \n) + - j drdr'n (r) ^ n (r') + F^e(/3 | n). (21) 

2, J I r r' I 

With these definitions ( !T4|) becomes (with an overall factor of (3 included to make the 
equations dimensionless) 

/3/i(°)(r,/3 I n) = /3;.o(r) - (3e' [ dv'n (r') - /^ ^^J^f '^^ . (22) 

J |r — r| on (r) /? 

Again, this relationship does not hold for all n. Rather, it is an equation to be solved for 
n = riQ. Its solution can be addressed in two ways. The Kohn-Sham approach is suggested 
by combining Eqs. f lT^ . flTK]) . and (^T^ in a form equivalent to 

5FW(/3 I n 



^ = /i - VKs{r,/3 I n) , (23) 



6n (r) 

where the KS potential vks is identified as 

VKs{r,/3 I n) = ^;o(r) + /" dr'n (r') + ^:^££i^ . (24) 

J |r — r I dn{r) i3 

Equation fl2^ has the form of a non-interacting Fermi system in the external potential 
Vks- The density and all other thermodynamic properties then are determined from the 
self-consistent single-particle Hamiltonian eigenvalue problem 

+ ""^^^"""^ ' " " ^ ' ''^^ ^""^'^ ■ 

Note that, in this approach, explicit construction of the functionals F^^\(3 \ n) and jj^'^\r,(3 \ 
n) is avoided, at the cost of introducing the KS orbital manifold, an issue already mentioned. 
The alternative is to construct these functionals, then to solve (l22l) directly as an equation for 
n, and use its solution in ( l2Ti) to determine the free energy. This second method is OFDFT; 
it avoids construction of the Kohn-Sham eigenfunctions and eigenvalues. Of course, when 
properly implemented the two approaches are equivalent. 

In the context of OFDFT, the remainder of this discussion is restricted to determination 
of F*^°^(/3 I n) and [5^^^\r,[5 \n). K suitable determination of the corresponding exchange- 
correlation functionals is a separate problem not considered here. 
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III. FUNCTIONALS FOR THE NON-INTERACTING SYSTEM 



A. Functionals of the density 

Without Coulomb interactions, the Hamiltonian of Eq. (jl]) is a sum of single-particle 
Hamiltonians. Consequently, the evaluation of the grand potential ([7]) can be reduced to a 
single-particle problem. For completeness, the analysis is sketched in Appendix |Al with the 
result 

I ^) = -p~^Tr In + e'^^^-^^^))^ , /i(q) = /i - t;(q) . (26) 

The trace is now taken over the single-particle Hilbert space (spin is not considered here for 
simplicity). The density is obtained by ([H]) (see also Appendix lAl) 

n(r,/3 I /.) = -^^^l^^^fM =rr5(r-q)p„ (27) 

OyU (r) /3 

where the one-body reduced equilibrium statistical operator (hereafter referred to as the 
"Fermi operator") is 

% = (^e-^(^-^(^» + ' ■ ^^^^ 
The central observation at this point is that ( 127|) defines the one-to-one relationship of /i(q) 
to an associated density 

/^(q)=/iW(q,/3|n). (29) 

A special case is the relationship for a uniform external potential for which /i*^°^(q, (3 \ n) ^ 
fi^e\(3,n), which is the familiar relationship of chemical potential to density for the ideal 
Fermi gas. The functional /i^°-'(q, /3 | ■) is universal, i.e. it is the same for all external 
potentials. Use of (129|) in (128|) expresses the Fermi operator as a universal (operator- valued) 
functional of the density as well 

= p,(/3 \n)=( e/^(^-^'"'(^i'/^l")) + iV • (30) 



The free energy for the non-interacting system is obtained from il^'^^((3 \ p) by the Leg- 
endre transformation, Eq. ffTOl) . as 



/3F(°)(/3 I n) = -Tr In 1 + e'^^^-^^'i" - /3p (q) pe . (31) 
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where use has been made of ( !27|) for the density. This result can be expressed in terms of 
the Fermi operator by ehminating (q) using fl28|) to get 



Tr 



[1 - Pe) In (1 - Pe) + Pe In Pe + Pe 

2m 



(32) 



Finally, entropy and internal energy functionals can be identified from the definition 



^(0) ^ ^(0) _ 2^5'(o) 



to get 



n] 



U^'^\P I n) 
-ksTr [(1 - Pe) In (1 - pe) + pe Inpe 



(33) 

(34) 
(35) 



Since the spectrum of pe is real and in the interval (0, 1), the energy and entropy are positive. 
Observe that all functional dependence upon the density occurs through the Fermi operator 
in the form of Eq. fl30l) . 

This development delineates the essential task, which is to expose these definitions as 
explicit functionals of the density. Consider, for example, the determination of p^^\q_, /3 \ n) 
from inversion of the density equation (1271) . Evaluating the trace in coordinate representation 
gives 

(36) 



r 



This expression is deceptively simple, since extraction of p'-'^-* is complicated by the non- 
commutation of p*^''''(q, /3 I n) with the kinetic energy. A common approximation is the 
Thomas- Fermi form, which is obtained by replacing the operator p^^-* (q, /3 | n) by the c- 
number p*^°^(r, /3 | n), so that the Fermi operator becomes diagonal in momentum represen- 
tation 



niv) — )■ h' 



dp ( e 



,/3{i^-/.(0)(r,/31n)) 



+ 1 



(37) 



In this case, the inversion is the same as that for the uniform ideal Fermi gas noted above 
p(°)(r,/3 I n) ^pi°^(/3,n(r)). This inversion is exact in the uniform limit, but conditions for 
its validity otherwise are not known. In the following discussion, several exact properties of 
the density functionals are exposed to provide guidance and constraints on approximations 
for their explicit forms. 
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B. Functionals of the statistical operator 



To put the equilibrium functionals of the last section in context, it is useful to define re- 
lated functionals of trace class operators p which are Hermitian, positive, and have spectrum 
on (0, 1) but otherwise are arbitrary. For them, define 



\p) = Tr 



^2 

:i-p)ln(l-p)+plnp + /3|-p 



(38) 



[p] = Tr£p , (39) 

5(0) [p] = -fcsTr [(1 - p) In (1 - p) + plnp] . (40) 

Clearly these functionals become the density functionals of the preceding section when eval- 
uated at the Fermi operator, p = pg. More generally, J^^^\[3 | p) is a convex functional of p 
while iS'-O) [p] is positive and concave. Furthermore, the minimum of J-'^^\(3 \ p) occurs for 
p = Pe, among all p with the same density 

n{r) = Tr(5(q — r)p . (41) 

To see this evaluate 

6 (^J^^°\f3 \P) + J dr\ (r) (n(r) - Tr6{q - r)p)^ = , 
where A (r) is a Lagrange multiplier function. The solution is 



.2 N N -1 



(42) 



Pm=|eV ^ P""- v43) 

The last equality follows since A (q) must be determined from Eq. fl41l) , which is the same 
as ( |27|) . In summary, 

I n) = I Pe) < I P)- (44) 

This result will be used below to establish upper bounds on the free energy density functional. 



IV. SCALING RELATIONS, BOUNDS, AND INEQUALITIES 
A. Scaling from dimensional analysis 

In this section scaling properties of the functionals are obtained from dimensional analysis. 
An alternative constructive derivation of scaling properties is given in Appendix [Bl It is 
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closer in formulation to the uniform scaling widely used in ground-state DFT 15(] but does 
not expose the simplicity and profound connection with thermodynamics that comes from 
the dimensional approach. 

All the density dependence of the density functionals occurs through the dependence of 
the Fermi operator on /i'^°)(r,/3 | n) or, equivalently, the dimensionless local "activity" 



z/(°)(r,/3 I n) =/3/i^"^(r,/3 | n) . 



(0)/ 



(45) 



Define the dimensionless momentum and position operators 



(46) 



where A is the thermal de Broglie wavelength, 



m 



2\ 1/2 



(47) 



Then Eq. f p6|) becomes 



n*(r*) = A'^n(r) = ( r 



(48) 



This shows that z/(r) is a functional of the dimensionless density n*(r*) alone 

z/(r) = z/(°)(r,/3 | n) = u'^°\r* \ n*). (49) 

The immediate consequence of (H9l) is that the Fermi operator and all functionals in dimen- 
sionless form are also functionals of n* only, to wit 

f3F^'^\P\n) = F*[n*], (3U^'^\p \ n) = U*[n*l -^S^''\p \ n) = S*[n*] . (50) 



Consider now the scaling transformation 



r — )-7r, n{r)dr — )■ 'y^n{'jv)dr, V ^ 'j "^V 



where V is the volume. This transformation preserves the total number of particles 



n{'r)dr = N = / 'j^n{'jr)dv . 

V J-y-'-iV 



The dimensionless density then transforms according to 



n*{r*) = X'n{r) 7^A'n(7r) = (7A)^n(7Ar*) . 



(51) 



(52) 



(53) 
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This is the expected result, namely that the length scale A used to define the dimensionless 
density itself has been scaled by 7. However, this scaling can be canceled by a suitable 
scaling of the temperature in the de Broglie wavelength, to wit 

7A(7-'/3) =A(/3) . (54) 

Thus, a combined scaling transformation of lengths and of the temperature leaves the dimen- 
sionless density invariant, and hence all the functionals of f l50|) . unchanged. To be explicit, 
define the transformed density by 

n^{r) = 7^n(7r) . (55) 
Then from fISU]) and the invariance of n*, the desired scaling properties follow: 

iP\n)= 7-2f(°) (7-2/3 \n^), f/(°) (/3 | n) = 7-2f/(o) (7-2/3 | n^) (56) 
5(°)(/3 I = 5(°H7-2/3 I n^). (57) 

B. Bounds and inequalities 

In this section specific upper and lower bounds are described. In addition, the scaling 
laws are used to obtain inequalities for the various functionals for different values of their 
arguments (temperature and density). 

1. Lower bound 



It is well-known [18|-[20|] that the non-interacting kinetic energy density functional at zero 



temperature is bounded from below by the von Weizsacker functional 



U^n] ^ — / dr^-^ = ^ / |Vvp(r)|^ . (58) 
8m J n[r) 2m J 

The last equality is obtained by introducing the square root of the density (a positive 
function) 

n(r) = (^(r)2 . (59) 

We now show that the von Weizsacker functional is a lower bound for the kinetic energy 
density functional at finite temperature as well. 
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Since pe is a Hermitian, positive operator it has a square root: 
Then U^^"^ {(3 \ n) from Eq. flM|) can be written 



[/(o)(/3 I n) = —Tr 



2 _ 

2m 



drdr' 



V (r| ^ |r') 



(60) 



(61) 



Note the similarity of this exact result to the approximate von Weizsacker result fl58|) (the 
relationship of il) to the square root of the density is '^{j) = \J (r| t/^^ |r)). To extract the von 
Weizsacker contribution, define 



X(r,r') = i^-^{Y) (r|^|r') 



(62) 



Then 



U^^)(l3\n) = — [ drlVMf [ dv" \X (v,r")f + — [ drlMf [ rfr"|VX(r,r") 
2m J J 2m J J 



+ 



2m 



drip{r)Wipir)V J dr" \X (r, r" 



But 



so 



dr"|Xfr,r")l' 



n r 



dr" (r| \r") {r"\ ■?/' |r) = 1 



(63) 



(64) 



f/W(/3 



n] 



J dr|V<^(r)|V|^ J dr\ip{r)\^ J rfr"|VX(r,r") 



2m _ 

> U,w M (65) 

This is the desired result, showing that U^w is a lower bound for U^^\ The analogue 
in ground-state DFT is the so-called Pauli-term decomposition |32|-|35| of the KS kinetic 

energy 

rs = %w + re , re>0. (66) 

Eq. ([nS]) also provides a bound for the free energy functional, namely 



F{(3 I n) > U,w M - TS{f3 \ n) . 



(67) 



In obtaining fl65|) . no use has been made of the integral over r, so provided that no 
traceless forms are added to the energy densities, the bound holds for them too: 



F{(3 \n) = j drf{r,(3\ n), U{(3 \ n) = j dru{r,l3 \ n), S{(3 \ n) = j drs{r,(3\ n) 



(68) 
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j dv"\VX (r,r")|' > u,w{r I n) 
f{r, (3 \n)> Uyw{r \ n) - Ts{r, f3 \ n) . 



u{r,j3 I n) = u^wir \ n) + - — n[r 

2m 



(69) 
(70) 



2. Upper bound 

Upper bounds are provided by if a suitable choice for p can be found. In practice, it is 
common to use the local density approximation, which in this context leads to the Thomas- 
Fermi form for the density. The question arises as to what statistical operator corresponds 
to the Thomas- Fermi approximation. More specifically, what p has the Thomas- Fermi form 
for the density? The necessary condition is 

(r \ptf\ r) = h-^ j dp |^e^(^-^=(/3'"W)) + . (71) 

However, this condition determines only the diagonal matrix elements of ptf- To fix the sta- 
tistical operator definitively, consider its associated Wigner distribution utf (p, R), defined 

by 

(r| PTF W) = 1 d^ei^-'^UTF (P, R) , R = q = r - r' . (72) 

Specialization to q = and comparison with ( ITTl) determines utf {p, R) as 



riTF (P, R) = (^e^(i^-'^^(^'"(^))) + 1 j . (73) 
With the statistical operator determined, an upper bound is given by (jUj) with p = Ptf 

(/3 I n) < ^(°) (/3 I ptf) = U^'^ [ptf] - [ptf] ■ (74) 
The calculation oUA^'^^ptf] is straightforward 

druTFi(3,n{r)) , (75) 
where UTF{P,n{r)) is the finite temperature Thomas- Fermi energy density 

UtfW, n (r)) = 1 1 rfp|^ |^e'^(i4-Me(/3,n(r))) (7g) 

To complete the calculation of the upper bound it is necessary to evaluate 

5(0) [ptf] = -ks jdr{r\ [(1 - ptf) In (1 - Ptf) + Ptf In Ptf] \ r) . (77) 
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Although the density dependence of all matrix elements of pxp is explicit, the evaluation of 
non-linear functions of this operator is still a difficult problem and the explicit form for f l77|) 
has not yet been determined. 

In any event, the upper bound for the free energy density functional is 

\n)< j druTF{(3,n{r)) 

+ (3-^ J dr{r\ [(1 - ptf) In (1 - Ptf) + Ptf In Ptf] I r) . (78) 

3. Inequalities 

A class of inequalities can be obtained from Eq. (1441) and the scaling laws. Recall the 
definition of in (!55|) . and choose p = 'pe{'~f~'^l3 \ n^) in (jH]), 

{/3 I n,) = (/3 I MP I ^7)) < J"^'^ (/3 I Mr'P I n,)) , (79) 
and then use the definition of J-''-'^^ (/3 | p) to get 

(/3 I n^) - (7-2/3 I n^) < (7^ - 1) TS^^^ (7-^/3 | n^) . (80) 

Use of the scaling laws (15611571) to convert functions of ■j~'^(3 to corresponding functions of /3 
leads to 

(/3 I n^) - 7^^ {/3 \n) < (7^ - l) TS^^^ {/3 \ n) . (81) 

Since the entropy is non-negative, an inequality for constant temperature and different 
densities is 

{(3 I n^) < 7'f(°) (/3 I n) , 7 < 1 . (82) 

Use of the scaling laws once more to convert (/3 | n^) to 72f(°) (7^/? | n) gives the in- 
equality for constant density 

(7^/3 I n) < (/3 I n) , 7 < 1 . (83) 
Repetition of this analysis for the choice 

(7-^/3 I n,) = (7-^/3 I Pe(7-'/3 I n,)) < (7-^/3 | p,(/3 | n,)) , (84) 

leads to 

(7-2/3 I n^) - (/3 I n^) < (l - 7^) T^(°)(/3 I ^7) , (85) 
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and scaling gives the reverse inequalities for 7 > 1: 

72^(0) (/3 I n) < (/3 I n^) , 7 > 1 
In consequence, it also is true that 

{(3\n)< (72/3 I n) , 7 > 1 . 



(86) 



(87) 



Appendix [C] has further consideration of Eqs. (180!) and (185|) to reach corresponding in- 
equalities for the entropy and energy functionals. The results are summarized in the Table. 
These inequalities imply that the entropy and energy are monotonically increasing with T, 
while the free energy is monotonically decreasing with T, at constant n . 



Constant n 


Constant /3 


(/3 


n) 1 (72/3 


n), 7 1 1 


F(o) iP 1 


n,) 1 7'i^(°) (/3 1 n) , 7 1 1 




n) 1 ?7(o) (7^/3 


^) , 7^1 


?7(o) (/3 1 


n^) 1 7'f/^°^ (/3 h) , 7 1 1 


5(0) (/3 


n) 1 ^(o)(72/3 1 


n), 7^1 


^(o)(/3 


|nA)g^W(/3|n), 7|1 



Return to (IHTj) and use scaling to write it as 

(7^/3 \n)-F{P\n)< 7 ^V gW (/3 | n) . 



7^ 



Set 7 = 1 + (5 with 5 — )■ to get 



n) 



(88) 



(89) 



This is a thermodynamic identity, which confirms the fact that the functionals considered 
here are all derived from the grand potential fl26|) for a non-uniform system at equilibrium. 



V. DISCUSSION 



The analysis of density functionals relevant to finite temperature DFT has been placed in 
the context of equilibrium statistical mechanics for a non-uniform system. The free energy 
functional of f lTU]) is precisely the same as that introduced by Mermin, but obtained here 
as a thermodynamic change of variables via Legendre transformation. This recognition 
of the thermodynamic context for DFT is expected to play a more important role when 
the functionals are parameterized by temperature, than for familiar ground-state DFT. 
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The origin of the density dependence through inversion of the relationship of n(r) to /i(r) 
allows a formal representation of all functionals for constructive evaluation: the more general 
constrained search methods are replaced by an explicit, albeit very difficult, many-body 
evaluation. The analysis here is restricted to the non-interacting functionals for which that 
many-body problem can be solved. Nevertheless, the necessary inversion to find /i(r) = 
/i''°^(r, /3 I n) remains a difficult mathematical task. Hence, construction of approximate 
functionals is required in this case as well. The objective here has been to provide both 
context and constraints for that process. 

The non-interacting functionals of section IIIII are deceptively simple and still defy exact 
evaluation. However, they are in a form convenient for exposing some of their general 
properties. Linear coordinate scaling is an example. For it, simple dimensional analysis 
provides the exact results. The results are analogous to those obtained for ground-state 
DFT, except the temperature is an additional variable that must scale along with the density. 
This is an example of a general outcome. Many of the finite-temperature results conform 
to obvious expectations from the ground-state theory, but some are not predictable by such 
extrapolation. 

Application of these results includes determination of the volume dependence of the func- 
tionals (see, for example, Chihara and Yamagiwa [13|). Identification of the von Weizsacker 
lower bound is perhaps expected, but the proof also provides a representation for the exact 
positive correction that may be suitable for further analysis. Upper bounds are provided 
by noting that all functionals depend on the density through the Fermi operator, and that 
the free energy is the extremum of a convex functional of constrained single-particle reduced 
statistical operators T^^^/S \ p). As an application, the single-particle reduced statistical 
operator corresponding to the TF approximation was constructed and shown to yield the 
TF approximation for a class of functionals, but not for the entropy in particular. Thus 
TF does not result from some approximate single-particle statistical operator. This will be 
discussed further elsewhere. 

The inequalities of the last section are a consequence of the scaling laws and convexity 
of the I p). Since the scaling transformations involve both temperature and density, 

two classes of inequalities are described: those for functionals at different densities but the 
same constant temperature, and for the temperature dependence of functionals for constant 
density. Ultimately, it is expected that such inequalities have an instructive thermodynamic 
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interpretation (see for example (189|) ). 



VI. ACKNOWLEDGEMENT 

This work was supported under US DOE Grant Ide-sc0002139l We acknowledge informa- 
tive conversations with Stefano Pittalis. 

Appendix A: Determination of I fi) 

1. Derivation of Eq. (l26]) 

Development of the thermodynamics for an ideal gas occurs in every textbook on equi- 
librium quantum statistical mechanics. For completeness, the extension to non-uniform 
systems to derive ( 126|1 is described briefly here. First write Eq. ([7]) in the form 

oo 

Qi^){l3 I ^) = _/3-iln^Tr(^)e-^(^"-"(5-P')) (Al) 

Af=0 

with the single-particle operator a(qj, pj) given by 

«(qi, Pi) = 77^ - /^(q*) (A2) 
2m 

Since a(qj, pi) is a Hermitian operator, it has an associated complete basis set for the single- 
particle Hilbert space 

a(qi. Pi) I ai >= Ea, | a, > . (A3) 

Observe that the "a" are state labels while the "i" are particle indices. The N particle 
Hilbert space then is spanned by the anti-symmetrized product basis S \ ai, ..,aN >, and 

N / ^ \ 

a(qi, Pi)S I ai, .., Oat >= I ^ J 5 | ai, .., oat > . (A4) 

i=l \i=l / 

These total energy eigenvalues can be specified in an equivalent way in terms of the occu- 
pation numbers for a given quantum state 

N 

= ^£ana ■ (A5) 
i=l a 

The sum on the right side is over all single-particle states, and is the number of times 
the quantum number a occurs in the state | ai, .., >, restricted by ^c^n^ = N. 
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Since the states are anti-symmetrized, there is no meaning to assignment of which par- 
ticles are in the state with quantum number a, only the total number of particles in that 
state ria- Consequently, the states can be labeled by {ria} instead of {aj} 

S I Oat >^ — 5-1 {no} > (A6) 

In this occupation number representation, the effects of symmetrization are accounted for 
since one set of {11^} represents all the states generated by permutations. The trace for the 

grand potential can be carried out in this occupation number representation as follows: 

00 

I ^) = -p-^\nY^J2 < K} I e-^(^-i"(^-P')) I {n^} > 

N=0 {n„} 

00 / 00 / \ 

= 1^ E n E = 1^ E n E ^"'"""^ E (AT) 

N=0 a Ua N=0 a Ua \ a / 

Primes over summations over ria indicate the restriction to Y2a = N , a. restriction that 
is removed in the last equality by the explicit Kronecker delta. The summation over is 
then performed first, giving 

Fermion occupation numbers are all = 0, 1 so 

I ^) = -l3-^Y^\n{l + e-^'-) = ^ In (l + e-'^'^^^^'P^)) | ai > 

a a 

= -/3-^Tr In (l + e-^^'^^''^'^) (A9) 
This is the result quoted in the last line of 



2. Proof of n(r) = -Sn'-^^ (/3 | /x) /6fi (r) |^ 

The functional derivatives entail use of some operator algebra. To illustrate this, the 
derivation of Eq. flTTI) from fl2Bl) is described. Begin with the functional derivative of 



djj (r) /3 (r 



r^Tr«T^ In 1 + e^'l^^'^^^V I . (AlO) 



To calculate the derivative of the logarithm a useful representation is 

•1 /.I 



ln(l+X)= [ dXdx\n{l + XX) = ! d\{l + \X)-^X , (All) 
Jo Jo 
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where 

The functional derivative is then 
6 



X = e 



/3U;7-M?,/3|n) 



(A12) 



6n{r) 



ln(l + X)=/ dX{l + XX)-'——{l + XXy' (A13) 
Jo on[rj 



In the present context, this expression occurs under a trace, so use of the cychc invariance 
of the trace gives 



5/x(r) 



p-^Tr [ dX{l + XXy 
Jo 



, 5X 

6n{r) 



r^Tr{l + X) 



.1 ^X 
5n{r) 



Next, the derivative of the exponential flA12p can be defined as 

6e 



(5/i(r) 



lim 



— e 



(5/i(r) 



Then writing 



an equation for is obtained: 



Integrating this equation gives 







1+ / rfx'e 
'o 



order 



and so 







The derivative is then 



6e 



(A14) 



(A15) 



(A16) 



(A17) 



(A18) 



+ order {SfiY . (A19) 



(A20) 
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With these results, and again using the cychc invariance of the trace, (1A10|) becomes 



5/i (r) 

which is Eq. (I27|) of the text. 
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l3-^Tr{eV^ ''^'V + l)-^5{v-v) = n{r) (A21) 



Appendix B: Scaling from unitary transformation 

Define the one-parameter family of operators 

U{z) = e^"(5-P+P-5) , (Bl) 

and the related operators 

X{z) = U{z)q_U-\z), Y{z) = U{z)pU-\z) . (B2) 
Then differentiation with respect to z 

^ = [(q-p + p-q) , q] U{z)-' = X{z) , (B3) 

^ = [(q-p + P-q) , P] U{z)-' = -Y{z) , (B4) 

gives (noting X(0) = q, 1^(0) = p) 

U{z)-^U-\z) = e"q, U{z)pU-\z) = e'^p, -oo < z < oo . (B5) 

One more definition is needed: 

= U{\n (7)) = ln(7){q P+P-q) _ ^Bg) 

This is the unitary operator that generates scale transformations 

U^A{^,p)U;' = A(7q,7^ip), < 7 < 00 . (B7) 

Now consider the density given by ( 127|) . insert unity U^^U^ = 1, and use cyclic invariance 
of the trace 

n(r) = Tr6{q - r)pe = TrU:^^U.f6{q - r)p = Tr6{'yq - r)pe^ , 
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or 

(r) = 7^n(7r) = Tr5(q - r)pe^ . (B8) 



Here is the transformed Fermi operator 



-1 



Pe-, = f/^Pet/-^ =\^e V- ' ' V + Ij (B9) 

Pe^ = f/^p,(/3 I = Pe(7-'/3 I n^) . (BIO) 

Uniqueness of the functional fj.^^\ci,l3 \ n) gives the identification 

/3/.(7q,/3 I n) = (7~2/3) p(q,7"'/3 I ^7) • (BH) 

Consequently, the statistical operator pe obeys the scaling relation 

Pe^(/3 I n) = Pe(7"'/5 I n^). (B12) 

Confined systems can be represented by including a wall potential in the external po- 
tential. For example, spherically symmetric confinement is represented v^, = Q {R — q) e^, 
where is the potential barrier for escape from the volume V = AttR^/S, with the barrier 
taken arbitrarily large for the hard- wall limit. Under the scale transformation, the wall 
potential becomes 

U^Q {R - q) U-^ = Q{R-^q) = e [j-^R - q) . (B13) 

Thus the volume transforms as — ?■ 7~^l^. 
Consider next the free energy functional (|32|) 



1 - Pe^) In (1 - pe^) + Pe^ lu Pe^ + Pl~'^7^Pe-i 



{P\n) = (3-'Tr^'^U;% [..] 

= 7-2^(0) (7-2/? I n^) (B14) 
The scaling for the energy density and entropy are found in a similar way 

f/(o)(/3 I n) = 7-2t/(°H7"'/3 I ^7), S^^'^/S \ n) = S^^^T"'/? I ^7) (B15) 
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Appendix C: Inequalities from scaling 

Equations (IHO!) and (185|) can be combined to give 

(1 - 7') T5W(7"'/3 I n,) < (7-2/3 I n^) - {/3 \ n,) < (l - 7^) TS^''\/3 \ n,) , 

(CI) 

and application of the scaling laws gives 

I n) < (/3 I n) - (72/3 | n)) < il^^TS(°)(7'/3 | n) . (C2) 

This implies 

^W(/3 I n) I 5(°)(7'/3 I n), 7|l (C3) 

or, equivalently, 

5(°)(/3|n)|5W(/3|n,), 7|l. (C4) 
Next, writing f l80p in terms of the energy and density 

f/(°) {13 I n^) - (7-2/3 I n^) < T (5(°)(/3 I ^7) - 5(°)(7"'/3 I %)) , (C5) 
or with scaling 

f/(°) (/3 I ra^) - 72f/(°) {I3\n) <T (S''^°^(/3 | ra^) - S'^^\l3 \ n)) . (C6) 
Then (EH]) and (El give 

[/(o) (/3 I n,) I 7't/^°^ (/3 I n) , 7 I 1 (C7) 
Expressing these results in terms of the same density by scaling gives the equivalent forms 

(7^/3 I n) g f/(°) (/3 I n) , 7 | 1. (C8) 
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